Cosmological Constraints from Type la Supernovae Peculiar Velocity Measurements 
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We detect the correlated peculiar velocities of nearby type la supernovae (SNe), while highlighting 
an error in some of the literature. We find <rg = 0.79 ± 0.22 from SNe, and examine the potential of 
this method to constrain cosmological parameters in the future. We demonstrate that a survey of 
300 low- 2 SNe (such as the nearby SNfactory) will underestimate the errors on w by -~35% if the 
coherent peculiar velocities are not included. 
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The first compelling evidence that the Universe is un- 
dergoing a period of accelerated expansion was provided 
by observations of Type la supernovae (SNe) d 0. 
The data from many current 0, 0, [B[ and near fu- 
ture [10, HQ surveys should eventually constrain the ef- 
fective dark energy equation of state to better than 10%. 

Density inhomogeneities cause the SNe to deviate from 
the Hubble flow, as gravitational instability leads to mat- 
ter flowing out of under-densities and into over-densities. 
These "peculiar velocities" (PVs) lead to an increased 
scatter in the Hubble diagram, of which several studies 
have been made [II El El El El, El El E3, El, El, HI 
l2lL I22I ] . When combining low and high redshift SNe in 
order to estimate the properties of the dark energy, the 
velocity contributions are usually modeled as a Gaussian 
noise term which is uncorrelated between different SNe. 
However, as recently emphasized [23], in the limit of low 
redshift z < 0.1 and large sample size, the correlations 
between SNe PVs contribute significantly to the overall 
error budget. In this letter we investigate the effect of in- 
corporating these correlations using the largest available 
low redshift compilation [17] . 



I. PECULIAR VELOCITY COVARIANCE 

The luminosity distance, di, to a SN at redshift z, is 
defined such that T = -rAr where T is the observed flux 

and C is the SN's intrinsic luminosity. Astronomers use 
magnitudes, which are related to the luminosity distance 
(in megaparsec) by 



fj, = m — M = 5 log 10 d,L + 25 



(1) 



where m and M are the apparent and absolute magni- 
tudes respectively. In the context of SNe, M is a "nui- 
sance parameter" which is degenerate with log(i/o) and 
can be marginalised over. For a Friedmann-Robertson- 
Walker Universe the predicted luminosity distance is 
given by 



d L (z) = (1 + 2!) 



H(z') 



(2) 



in speed of light units, where H is the Hubble parameter. 
In the limit of low redshift this reduces to dh ~ z/Hq. 

Assuming a flat Universe (as we shall do throughout 
this letter), the effect of peculiar velocities leads to a 



perturbation in the luminosity distance (<5di) given by 
MI2IIMI2EI27] 



5dL 



= r • 



(1 + ^) 2 
H(z) d L 



v - v _ 



(3) 



where r is the position of the SN, and vq and v are 
the peculiar velocites of the observer and SN repectively. 
Using the Cosmic Microwave Background dipole we can 
very accurately correct for vq. This demonstrates how 
a SNe survey that measures fj, and z can estimate the 
projected PV field. We now relate this to the cosmology. 

The projected velocity correlation function, £(ri,rj) = 
((v(i"i) ■ fi)(v(rj) ■ fj)), must be rotationally invariant, 
and therefore it can be decomposed into a parallel and 
perpendicular components [28l. |29|. [30j : 

£(ri,rj) = sin 6^ sin6j£±(r, zu z,)+cos 6^ cos 6>j£|| (r, z,-, Zj) 



where ry 



ri i-j , r = |ry |, cos 6i = r i • f y , and cos 9 
In linear theory, these are given by [H, [29|, [3; 



ill 



, ± =D\z l )D'(z J ) ^P(k)K 



1,-L 



(fcr) 



(4) 



where for an arbitrary variable x, Ku (x) = ja(x) — 2jl ^ , 
Kx(x) = ji(x)/x. D(z) is the growth function, and 
derivatives are with respect to conformal time. P(k) is 
the matt er p ower spectrum. This corrects the formulae 
used in [HHH, see the Appendix for details. 

The above estimate of £(ri,r_j) is based on linear the- 
ory. On scales smaller than about 10/i _1 Mpc nonlin- 
ear contributions become important. These are usually 
modeled as an uncorrelated term which is independent of 
redshift, often set to a v ~ 300 km/s. Comparison with 
numerical simulations [33[ has confirmed that this is an 
effective way of accounting for the nonlinearities. Other 
random errors that are usually considered are those from 
the lightcurve fitting (fA e rr), and intrinsic magnitude scat- 
ter (cr m ) found to be ~ 0.08 in the case of [l7| . It is just 
these three errors that are usually included in the analy- 
sis of SNe. 

In Fig. [1] we compare the covariance of bdhjdj, from 
peculiar velocities for a pair of SN 



C v (i,j)= 1 



Hdr 



(1 + ^ 
Hdr 



£(ri,rj) (5) 
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FIG. 1: The ratio of the covariance from peculiar velocities 
C v , compared to the random errors a, for a pair of supernovae, 
over a range of angular separations. In the upper panel we 
vary the separation on the sky, 0. In the lower panel we 
vary the redshift, with them both at the same z or with one 
supernovae fixed at z = 0.03 (dash-dot). 



to the standard uncorrelated random errors given by 



a« 2 = 



« + /w(ir)+ 1- 



(1 + *) S 
Hd L 



(6) 



with {n m ,CL b ,h,n s ,w,a s } ={0.3,0.05,0.7,0.96,-1,0.85} 
and {/x err , a m , a v } — {0.1,0.08,300}. We see that the 
PV covariance is comparable with the uncorrelated er- 
rors for low-z, and we note that the correlated errors 
become more significant the larger the dataset. 



II. CONSTRAINTS FROM CURRENT DATA 

A uncorrelated-error-only analysis of SNe allows one 
to constrain the cosmological parameters through the 
Hubble parameter in @. Namely tt m and w, for a flat 



Universe. However, by fitting for the PV covariance we 
probe the matter power spectrum and can therefore con- 
strain further parameters such as Ot. , Ho, n s , as, where 
these have their usual meaning. In the following analysis 
we also allow the SNe "nuisance" parameters M, a m , a v 
to vary. These are often set to fixed values, however 
marginalizing over them allows a better estimate of the 
uncertainty in the other parameters. The weighted in- 
tegration of the matter power spectrum in Eq. (|4|) has 
similar weights to those used in evaluating a$. There- 
fore, the higher <7g, the stronger the correlations between 
the projected velocities, obeying C v oc a\ for the other 
parameters fixed. 

We analyse nearby supernovae (z < 0.12) from [T3| 
who find improved luminosity distances to 133 super- 
novae from a multicolour light curve method. Follow- 
ing we exclude 9 supernovae from the set. These 
are supernovae that are unsuitable due to bad lightcurve 
fits, those who have their first observation more than 20 
days after maximum light, are hosted in galaxies with 
excessive extinction (A v > 2.0 mag) and one outlier 
(SN1999e). This leaves 124 supernovae z € [0.0023, 0.12], 
which have an average separation of f = 108 h~ 1 Mpc, 
a mean redshift z = 0.024, and herein we refer to this 
dataset as our "low-z" SNe. 

The likelihood is given by 

£ oc isr 1 /^-^" 1 * 

where — C v (i,j) + a(i) 2 Si 3 and x — [d°^ s — 

d L (z)]/d L (z) with df s given by Eq. {T}, and d L (z) by 
Eq. ([2]). The matter power spectrum can be evaluated ei- 
ther numerically (e.g. CAMB [34]) or using analytical ap- 
proximations [351 ] . We assume a flat Universe with a cos- 
mological constant (w = —1), a Big Bang Nucleosynthe- 
sis (BBN) prior Q b h 2 ~ AA(0.0214, 0.002) IH, and a Hub- 
ble Space Telescope (HST) prior h ~ A/"(0.72, 0.08) 
These two priors remove models that are wildly at odds 
with standard cosmological probes, but do not unduly 
bias results towards standard cosmology. The likelihood 
has almost negligible dependence on n s , and to keep it 
in a range consistent with CMB and large scale structure 
estimates we give it a uniform prior n £ [—0.9, 1.1]. The 
other parameters are all given broad uniform priors. We 
use the standard Markov Chain Monte Carlo (MCMC) 
method to generate sam ples from the posterior distribu- 
tion of the parameters [381 ]. Convergence was checked 
using multiple chains with different starting positions, 
and also the R — 1 statistic [39| . We also checked that 
the estimated posterior distributions reduced to the prior 
distributions when no data was used [H. All the anal- 
ysis was checked with two completely independent codes 
and MCMC chains. 

The low-z results are given in row A of Table HI As a 
non-zero trg is needed for the velocity correlations, these 
results indicate the correlations are detected at the 3.6a 
level. We also perform the MCMC runs without includ- 
ing the PV covariance matrix C v , and we find —2 ln£ max 
(where C is the likelihood) increases by 19.3. As the 
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likelihood no longer depends on {erg, fif,, h, n s } we have 
effectively removed four parameters. 

When estimating the cosmological parameters from 
SNe, a low redshift cut is usually imposed, to reduce the 
effects of the PVs. For example in [40(, their SNe dataset 
have 192 SNe with z E [0.016, 1.76] and a mean redshift 
of z = 0.48. Herein we refer to this dataset as our "high- 
s'' SNe. Our MCMC results for just this data (without 
including PVs) are shown in row B of Table HI Although 
we marginalize over the SNe parameters {M, a v , a s }, our 
constraints on U, m are still in excellent agreement with 
those obtained in [4(3]. 

We now combine the low-z and the high-z datasets, to 
make an "all-z" dataset. We used the overlapping SNe in 
the two datasets to estimate a small normalizing offset 
to the magnitudes from the latter data set, (the extra 
magnitude error is negligibly small) . The same procedure 
was used by in constructing their dataset. After 
eliminating duplicated SNe, our combined all-z dataset 
has 271 SNe with z £ [0.0023,1.76], and z = 0.35. Our 
all-z results are given in row C of TableQ] The constraints 
on <7g broaden slightly due to a mild degeneracy between 
us and Q m which is broken by the addition of the higher 
z SNe, but pushes as to the region of higher uncertainty. 
The increase in the minimum of —2 In/!, when as is set 
to zero, is now 16.4. 

One the main aims of cosmology is to test the cosmo- 
logical constant prediction w = — 1. Due to a degeneracy 
between w and f2 m , it is necessary to combine the SNe 
with another data source. Here we use the WMAP CMB 
data [4l|. We continue to assume a flat Universe, and 
we now allow w to be a free parameter. The results are 
given m rows D to H of Table Hand in Fig. [2] As can be 
seen by comparing rows D and E, if a redshift cutoff of 
z > 0.016 is made, then current data has a systematic er- 
ror of Sw = 0.02 when PVs are not included. This is sev- 
eral times smaller than the statistical error of Sw = 0.08. 
However, comparing rows F and G, shows that if no red- 
shift cutoff is made, neglecting the correlated PVs results 
in a systematic error of Sw = 0.07 which is about as large 
as the statistical error. 



III. FORECASTS 

We now consider the relevance of peculiar velocities to 
future supernovae surveys, using a Fisher matrix analy- 
sis. The Fisher (information) matrix probes the ability 
of an experiment to constrain parameters, by looking at 

the dependence of the likelihood: F a p = — ^ ® p \ 

F aP = d,„ C- 1 ^ +iTr(C- 1 C )Q C- 1 C, /3 ) (7) 

Often the second term is ignored, however we find that 
this approximation is no longer valid when including PVs. 
Unmarginalised la errors on parameter p a are given by 
>V (l/F Qa ), and the equivalent marginalised errors by 
y/ ({F _1 } QCt ). Generally, the inverted Fisher matrix F _1 



provides the expected covariance matrix for the parame- 
ters. 

We consider two future supernoave experiements that 
aim to detect high and low-redshift supernoave respec- 
tively: the Supernova/ Acceleration Probe (SNAP, [9|]), 
and the Nearby Supernova Factory (SNfactory, @). The 
baseline SNfactory program is to obtain spectrophoto- 
metric lightcurves for SNe in the redshift range 0.03 < 
z < 0.08, with the assumption that these SNe are far 
enough away for correlated peculiar velocities to not con- 
tribute significantly to the error budget. However, in 
Fig [1] we can see the contribution of correlated PVs is 
not irrelevant for z ~ 0.03, and as the number of SNe in- 
creases the overall uncertainty from random errors (such 
as intrinsic magnitude scattter and instrumental noise) 
gets beaten down by a 1/y/N factor, while the correlated 
errors from coherent peculiar velocities do not. Thus at 
any redshift, these peculiar velocity errors will begin to 
dominate for some large number of SNe. We investigate 
the situation for the SNfactory by considering 300 SNe 
randomly distributed over a rectangular area of 10,000 
sq degrees, and redshifts 0.03 < z < 0.08. We also in- 
clude high-z from a SNAP-like survey, and model this as 
2000 SNe randomly distributed over 10 sq degrees with 
0.2 < z < 1.7, and we assume PVs are irrelevant for these 
SNe. We take our fiducial model to be a flat ACDM cos- 
mology with fl m =0.3, fit — 0.05, h = 0.7, n s = 0.96, 
w = —1, ag = 0.85, and nuisance parameters a v = 300 
km/s, <r m =0.1, /J, err = 0.1. We further marginalise over 
M. 

In Fig [3] we compare the marginalised la contours ob- 
tained when the coherent PVs are ingnored and included, 
for our hypothetical SNAP and SNfactory surveys. We 
see that, even for a cut of z > 0.03 the error bars on O m 
and w will be considerably underestimated if the peculiar 
velocites are ignored, and in particular the marginalised 
error on w increases by 35% from 0.062 to 0.084, for 
the SNe alone. We therefore conclude that it is essen- 
tial to include a full covariance matrix analysis at these 
redshifts, to avoid significantly underestimating the er- 
rors. That the error bars increase rather than decrease 
indicates that the extra information available in the pecu- 
liar velocities is outweighed by the extra parameter space 
{Clb,h,n s ,a s }. 

In light of the extra effort required for the low-z SNe 
we consider excluding them. However in Fig [3] we see 
that the contours increase significantly when low-z SNe 
are not available (marginalised error on w increase to 
0.12), even when CMB data is included. Therefore we 
conclude that low-z SNe play a vital role in constraining 
the cosmological parameters, but a full covariance matrix 
analysis must be done. 

We have shown that a low redshift cut of 0.03 is futile 
in avoiding peculiar velocites when we have 100's of SNe. 
We therefore consider that no lower bound on redshift is 
needed at all, and we extend our study to include SNe 
with z < 0.03. We consider a survey of 1000 SNe ran- 
domly distributed over 10,000 sq degrees and redshifts 
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TABLE I: Mean and 68% confidence limits on the cosmological parameters, 
and a v ,a m , using different combinations of the WMAP and SNe datasets, with 
and without including the peculiar velocity covariance matrix (PV). See text for 
discussion. 




FIG. 2: 95% confidence limits on Q. m and 
w, for WMAP only (thin black), WMAP 
with high-2 SNe (thick red), and WMAP 
with all-2 SNe including PVs (dotted blue) 
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FIG. 3: The marginalised la contours for Q, m and w from 
a SNAP-like high-2 SNe survey in a flat ACDM cosmology 
(short dashed lines). We also consider including 300 low- 
2 SNe from a SNfactory-like survey, while ignoring peculiar 
velocities (solid lines) and including them (long dashed lines). 
Smaller red contours include cosmic variance limited CMB, up 
to I = 2000. 
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FIG. 4: The marginalised la contours for Q. m and erg, from 
1000 low-2 SNe in a ACDM cosmology, with 0.03 < 2 < 0.08 
(black solid line) and 2 < 0.08 (red dashed line), combined 
with a SNAP-like high-2 SNe survey. 



z < 0.08, and compare it to a similar survey of 1000 
SNe, but with 0.03 < z < 0.08. We combine them with a 
high-z SNAP-like survey as above. In Fig 0] we see that 
by including very low-z SNe (z < 0.03) we have gained 



free information about ag - 'free' because the error bars 
on tt m (and w) stay the same. However, a lower bound on 
z may be useful in reducing the effects of non-linearities. 

IV. DISCUSSION 

We have found that when no redshift cutoff is imposed, 
the current SNe data detects correlations in PVs at the 
3.6cr level. We also found that the current constraints on 
w are not sensitive to the PVs if SNe with z < 0.016 are 
excluded. However, we have shown that future large data 
sets will be sensitive to the PVs even if a lower bound of 
z > 0.03 is used. 

An alternative approach of accounting for the corre- 
lated peculiar velocities is to estimate the underlying den- 
sity field from galaxy redshift surveys and then use this 
to try and remove peculiar velocity at each SNe [2l[ . As 
this method has different systematics to the statistical 
modeling method we have investigated, we believe it will 
be a useful cross-check to apply both methods separately 
to future SNe data sets and compare the results. 

Appendix 

For perturbations, p(l +S), and peculiar velocity v we 
find the perturbed Fourier space continuity equation, in 
co-moving coordinates, is c) k — ik • Vk — where prime 
indicates differentiation with respect to conformal time, 
77. Using the linear approximation S^rj) = £>(?7)<5k, and 
assuming that the Universe has no vorticity (V x v = 0) 

leads to Vk = — i-D'^t-k, where D{z) is the growth func- 
tion. Therefore, the Fourier space correlation between 
the ith component of the velocity field at time tja and 
the jth component of the velocity field at time rjs is given 
by 

( V LX;) = D'( VA )D'( VB )(2nf5(k A -k B )P(k A )^, 

where (<5k A <5k B ) = (27r) 3 c>(kA — kB)-P(fcyi), and the sub- 
scripts denote the quantity at time tja or r\B- This cor- 
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rects equation (Bl) of 1311 . and our resulting Eq. ([5]) for 
C v corrects those of [32|, among others. 
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